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ABSTRACT 

Context. Field stars are not always single stars, but can often be found in bound double systems. Since binary frequencies in the birth 
places of stars, young embedded clusters, are sometimes even higher than on average the question arises of how binary stars form in 
young dense star clusters and how their properties evolve to those observed in the field population. 

Aims. We assess, the influence of stellar dynamical interactions on the primordial binary population in young dense cluster environ- 
ments. 

Methods. We perform numerical N-body simulations of the Orion Nebula Cluster like star cluster models including primordial binary 
populations using the simulation code nbody6++' 

Results. We find two remarkable results that have yet not been reported: The first is that the evolution of the binary frequency 
in young dense star clusters is independent predictably of its initial value. The time evolution of the normalized number of binary 
systems has a fundamental shape. The second main result is that the mass of the primary star is of vital importance to the evolution 
of the binary. The more massive a primary star, the lower the probability that the binary is destroyed by gravitational interactions. 
This results in a higher binary frequency for stars more massive than 2M Q compared to the binary frequency of lower mass stars. The 
observed increase in the binary frequency with primary mass is therefore most likely not due to differences in the formation process 
but can be entirely explained as a dynamical effect. 

Conclusions. Our results allow us to draw conclusions about the past and the future number of binary systems in young dense star 
clusters and demonstrate that the present field stellar population has been influenced significantly by its natal environments. 

Key words, binaries: general, galaxies: clusters: general, galaxies: clusters: individual: ONC, methods: numerical 



1. Introduction 

Today the formation process of isolated single stars seems to be 
understood in principle. One remaining uncertainty, at least to 
some extent, concerns the formation of binary stars. There are 
several ways to form binary stars - the fragmentation of a col- 
lapsing cloud or core or the dynamical capture of a companion. 
However, determining the dominating process of binary forma- 
tion requires knowledge of the parameters of the primordial bi- 
nary population, i.e. the population as it appears just after star 
formation has finished. 

Several surveys have been performed to determine the 
present-day parameters of binary populations in different kinds 
of environments. Observations of the field binary population 
have shown that the binary frequency is quite high, 60% for 
G-type dwarf stars (|Duquennoy & Mayor|1991[) and 35 - 42% 



for M-type dwarfs ( Fischer & Marcy 1992; Reid & Gizis|l 997). 
Both |Duquennoy & Mayor| ( |1991| l and |Fischer & Marcy| ( |1992| l 
determined the mass ratio and separation/period distributions of 
their samples. They found that the mass ratio distribution for G- 
type dwarfs increases as the mass ratios decreases while the mass 
ratio distribution for M-dwarfs is flat. The separation/period 
distributions for G-type and M-type dwarfs can be fitted by a 
Gaussian distribution in logarithmic space where both peak at 



nearly the same separation (~ 30AU for G-type and 9 - 30AU 
for M-type dwarfs) with the FWHM of both distributions span- 
ning over three orders of magnitude (see Duquennoy & Mayor 
[19911 Fig. 7). 

Another feature of binary populations is their dependence on 
the primary mass: observations indicate that the probability of 
finding a massive star in a binary system is higher than for low- 
mass stars. For OB stars, a binary frequency of ~ 70% is found 
( jPreib isch et al. 20 02|, whereas for G-d warfs a binary frequency 
of ~ 60% (Duquennoy & Mayor 1991 ) and for M-dwarfs a value 
as low as ~ 31% (March al et al.|200 3 ) has been observed. 

An understanding of the presently observed properties of 
evolved binary populations requires knowledge of the primor- 
dial binary population. To achieve this one has to look at sites 
where star formation has just happened, such as young star clus- 
ters. However, observations of binary populations in star clus- 
ters do not provide such a clear picture as the observations of 
the field binary population. While in dense star clusters such as 
the Orion Nebula Cluster (ONC) or NGC 2024 the obse rved bi- 
nary frequency is similar to the field binary frequency (Kohler 
|et al.|2006)|Beck et al.|20 03), sparse clusters such as the Taurus- 
Auriga association have much higher binary frequencies (see 
Duchene 1999 and references therein). This leads to the ques- 
tion of what defines the binary frequency in the different environ- 
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merits. Possible explanations are differences in the temperatures 
of the parent clouds ( Durisen & Sterzi kp"994 1 or the dynamical 
evolution within the cluster (e.g. Kroupa et al.|19 99). The latter 
possibility will be investigated in this study. 

According to current knowledge, stars do not form in iso- 
lation, but in a clustered environment consisting of several tens 
to thousands of stars. Most of them possibly formed in clusters 
such as the ONC or even denser environments ( |Lada & Lada| 
|2003||Pfalzner|2009| >. The stellar density in ONC-like clusters is 
sufficiently high to potentially affect the evolution of the binary 
population via gravitational interactions. 

The gravitational interactions that alter binaries in a star clus- 
ter are mostly three- and four-body interactions between binary 
and single stars. Heggie ( j 1975a|b) > and |Hills| ( | 1 975 1 > showed that 
three-body interactions are very effective in altering the param- 
eters of binary systems. Both authors found that loosely bound 
binaries (also called soft) tend to become softer (meaning that 
their internal energy increases resulting in the widening of the 
binaries) owing to three-body interactions, while firmly bound 
(hard) binaries tend to become harder (Heggie-Hills law). 

Four-body interactions (binary-binary encounters) become 
important in clusters that consist of a high number of binary 
stars (see for example Fig. 1 in |Leigh & "S ills 2010). In addi- 
tion to three-body interactions, four-body interactions can form 
two "new" binary systems by exchanging the binary compo- 
nents, can destroy one or even both binary systems, or can lead 
to the formation of bound hierarchies of three or more stars. 



NBODY6 ++ (|Aarsefh]|2003 Spurzem||1999|), the simulation code 



used for this study, uses several techniques such as the regular- 
ization of three -body interactions ([Aarseth & Zare|1974"| or the 
chain method (Mi kkola & Aarseth||1993 i to compute the evo- 
lution of these interactions correctly. However, in this study we 
only diagnose the evolution of the binary population. 

Numerical investigations of the dynamical evolution of bi- 
naries in star clusters started in the 70's with the development 
of the first N-body codes (e.g. Aarseth][l971 1. The early works 



concentrated on the formation of binaries in star clusters with 
initially pure single stellar systems. Only in recent years sim- 
ulations have been performed to investigate the evolution of a 
primordial binary population. For example, |Kroupa et al.| (l999 ) 
simulated the evolution of ONC-like star cluster models where 
all stars initially were part of a binary system, i.e. a primordial 
binary frequency of 100%. They showed that the cluster evolu- 
tion destroys wide binaries quite efficiently leading to a reduc- 
tion in the initial binary fraction as the cluster evolves. 

Here we also want to study the evolution of the primordial bi- 
nary population of a young dense star cluster. In contrast to the 
previous studies, we investigate the effect of the initial binary 
frequency (0 - 100%) and a log-uniform semi-major axis dis- 
tribution similar to that observed in very young star forming re- 
gions (e.g. |Reipurth et al.|2007 Connelley et al.[ 2008). Because 
observations indicate that the binary frequency increases with 
stellar mass, we also focus on the dependence of the binary evo- 
lution on stellar mass. 

This paper is structured in the following way. In Sec. [2] we 
define parameters that describe a binary population. Sec. 131 cov- 
ers the setup of the simulations, including the generation of the 
cluster models and initialization of the primordial binary popu- 
lation. Finally in Sec. [4] we present the results obtained from the 
simulations. 



2. Definitions 

The most simple measure of the binary population of a star clus- 
ter is the number of binaries {Nb) within the cluster. Here, a pair 
of stars is called a binary system if the system is bound and the 
two stars are at the same time nearest neighbors. In this defini- 
tion, a star that is part of a multiple system can be part of several 
"binary" systems simultaneously since two stars do not have to 
be mutually nearest neighbors. Hence triple or other higher or- 
der multiples are not handled separately by our diagnostic but 
are counted as multiple binary systems. In this study, we mostly 
use the normalized number of binary systems Nb defined as 



N b (f) = 



N b (t) 

N b (oy 



(i) 



where Nb(0) denotes the initial and Ni,(t) the number of binary 
systems at time t within the cluster. The usual measure of binary 
populations in star clusters is the binary frequency b / 



b f (t) = 



N b (t) 



N s (t)+N b (ty 



(2) 



where N„(f) denotes the actual number of observed single stars 
and Nb(t) the actual number of binary systems. One major result 
of our work is that the normalized number of binary systems, 
Eq. [TJ is the appropriate measure of a binary population. 

3. Method and cluster setup 

Here the evolution of binary populations in young dense star 
clusters is studied using the parallized version nbody6 ++ of the 



high-accuracy N-body code nbody6 (Aarseth 2003 Spurzem 



1999[). The study pre sented can be understood as a follow up 
of |Pfalzner & O lczak (2007). The main difference to Pfal zner &| 
|Qlczak| ( |2007| l is that here the clusters contain a primordial bi- 
nary population, with varying initial binary frequencies. In both 
studies, the ONC has been chosen as a model cluster because it 
is one of the most well-studied star forming regions for which 
many parameters are well constrained (see for example Prosser 



et al. 1994; Hillenbrand 1997; Hillenbrand & Hartmann 1998 



in particular several surveys of its binary popula tion (|Petr et a 



1998 



Preibi sch et al^[T999j |K6hler et al.1[2006t |Reipurth et a 



2007 ) exist. This allows a detailed comparison of the simulated 
binary populations with the observed ones. 

Our simulations started with 4000 stellar systems, where 
stellar systems correspond to either single stars or binary sys- 
tems. This means that for a simulation with a primordial binary 
frequency of 100% the dynamical evolution of 8000 stars is cal- 
culated. The masses M of the 4000 stellar systems were sampled 



from the initial mass function given by Kroupa ( 2002 ) 



f(M) 



M~ 13 
AT 2 - 3 
AT 23 



if 0.08 < M [M ] < 0.5, 
if 0.5 < M[M ] < 1.0, 
if 1.0 < M [M Q ] < oo. 



(3) 



The minimum and maximum masses of the stellar systems have 
been chosen as M m [ n = 0.08M o , corresponding to the hydrogen- 
burning limit, and M max = 5OM the mass of the most massive 
stellar system in the ONC - 6 l C Ori. 

Observations of the stellar content in the ONC have shown 
that the stellar radial density profile of the cluster can be approx- 
imated by an isothermal sphere (i.e.pn resent oc 1 /r 2 , see e.g. Jones 
& Walker| ( [T988| ; |Hillenbrand| ( |1997f ) with the exception of the 
cluster core CR C ore ~ 0.2pc), which has a much flatter density 
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profile (ppresent.core <* f " U5 , Scally et al. C 

this, we adopt an initial density profile of 



|2005 i). To account for 
the form 



Po(r) = 



-,-2.3 



A) (r/Rcoie) 

Po (r/R cole r 20 




r e (Q,R C0K ] 
re(R cole ,R] 
r € (R, oo] 



(4) 



where po = 3.1 x 10 pc , R CO ie = 0.2pc, and R - 2.9pc. Olczak 



|et al.| ( |2010[ ) demonstrated that this initial radial density profile 
evolves towards the currently observed one after about lMyr of 
dynamical evolution, the estimated age of the ONC ( |Hillenbrand| 
[l997] l. 

Observations of the ONC indicate that the cluster are mass- 



segregated (e.g. Hillenbrand 1997 1. Numerical simulations by 
Bon nell & Davies| ( |1998[ ) comparable to the cluster models stud- 
ied here demonstrated that this mass segregation cannot be due 
to dynamical evolution of the cluster and therefore must be pri- 
mordial. Following their suggestion, we initially place the most 
massive stellar system at the cluster centre and the three next 
most massive stars at random positions within a sphere of ra- 
dius R = 0.6Rhm centred on the cluster centre, where Rh m is the 
half-mass radius of the cluster. 

The velocities of the stellar systems are sampled from a 
Maxwellian velocity distribution in virial equilibrium. The re- 
sulting velocity dispersion of our mode l is lower than the obser 
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Fig. 1. Initial semi-major axis distribution for all simulated clus- 
ters. Crosses, stars, open and filled squares denote simulations 
with initially 30%, 50%, 75%, and 100% binaries, respectively. 



A flat binary period distribution has not only been observed 
in Sco OB II but also in other regions (e.g. 



vational value (<x = 4.3 + 0.5km s 

However, because of significant mass of the molecular cloud 
associated with the ONC that was neglected in our numerical 
model, the observed velocity dispersion is consistent without any 
assumption about a stellar system in dynamical equilibrium. In 
future investigations of the binary populations, we plan to in- 
clude the effect of the gas associated with the cluster. 

For the primordial binary population, we chose the model 
described in |Kouwenhoven et al.| ( |2007) . They compared obser- 
vations of the Scorpius OB II association with simulated obser- 
vations of their binary population model to eliminate observa- 
tional selection effects and biases that affect real observations 
of star clusters. Their analysis shows that the binary frequency 
in Scorpius OB II is above 70% with a confidence of 3cr. The 
closest agreement with observations is obtained for models with 



Connelley et al. 

2008; Reipurth et al. 2007) and has been used in several theo- 



Jones & Walker] [T988t . retical stu dies involving the evolut i on of primordial binary fre- 
quencies (|Heggie & Aarseth|1992 Ivanova et al.|2005 Hurley 



|et al.|2007D . 

In addition to a log-uniform shape of the initial period distri- 
bution, other distributions have been used in numerical studies 
of the dynamical evolution of binary populations. For example, 
Kroupa ( 1995a| l demonstrated that a dynamically evolved period 



a binary frequency of 100%. The model of Kouwenhoven et al 
(2007) was chosen here since it is a well-studied binary popula- 



tion, close to its initial conditions. Scorpius OB II is a relatively 
young (about 5 - 20Myr) and probably too sparse region to have 
undergone significant dynamical evolution. 

In contrast to the field population that exhibits a log-normal 



period distribution of G and M-dwarf binaries (Duquennoy & 



Mayor 1991| Reid & Gizis 1997 1, the younger population in Sco 



OB II has a flat distribution of the semi-major axis in logarithmic 
space ( Kouwenhoven et al.|2007 1 

f a (a)cc a - x o /iog 10fl (log 10 fl) = const., (5) 

where a is the semi-major axis of the binary system. 

A flat distribution of the semi-major axes in logarithmic 
space corresponds to a similar distribution of periods P with the 
relation 



1 



An 1 



logl0fl= 3 1Ogl ° P -3 1O °\GM r 



O"log l0fl - ^O-\ 0%m P, 



(6) 
(7) 



distribution that is initially rising for short periods and nearly 
constant for long periods matches the observed period distribu- 
tion of the field G-type dwarfs ( |Duquennoy & Mayo r|1991[ ). In 
contrast, this cannot be achieved with an initially log-uniform 
period distribution because of the much larger number of dynam- 
ically stable short-period binaries. However, physical processes 
that were not taken into account by Kroupa ( 1995a) could sig- 
nificantly lower the number of short period binaries. For exam- 
ple "the orbital decay of embedded binary stars" ( |Stahler|2 010) 
is capable of reshaping a log-uniform binary period distribution 
within less than 1 Myr towards a log-normal distribution (Hovel 
et al., in preparation). Consequently, as the effect of different ini- 
tial period distributions on the evolution of the binary population 
in star clusters is basically unknown we will explicitly focus on 
this topic in a future investigation. In this context the present 
study can be understood as preparatory work. 

Figure [TJ shows the initial semi-major axis distribution of 
our simulations. Although the semi-major axis distributions have 
been set up log-uniformly in the range 4.6 x 10~ 2 AU - 4.6 x 
10 4 AU, the distribution is only flat up to about 2 x 10 3 AU and 
declines afterwards to zero. The apparent lack of wide binary 
systems in the 2 x 10 3 - 4.6 x 10 4 AU range is due to the gener- 
ating algorithm. In the code, binary systems are formed in iso- 
lation and then placed in the cluster environment. The wider a 
binary, the higher the probability that (at least) one of its compo- 
nents has a closer neighbour among the other cluster members. 
According to our definitions (see Sec. |2]i, this chance neighbour 
is now the potential companion and most probably not bound. 

|Kouwenhoven et aT7| ( |2007[ ) also found the mass ratio distri- 
bution to have the form 



where G is the gravitational constant and Mj the total mass of a 
binary system. 



f q (q) «■ q y 



y q ~ -0.4, 



(8) 
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Fig. 2. Comparison of the stellar initial mass function (thin lines) 
and the initial binary-system mass function (thick lines). The 
open circles denote simulations that used random pairing to pair 
the binary components while the solid squares denotes simula- 
tions that used Eq. [8] 



where q = M sec /M pr j m is the mass-ratio of a binary with pri- 
mary mass M pl -i m and secondary mass M sec . Adopting this func- 
tional form in our models, the binary system masses were sam- 
pled from the |Kroupa] ( |2002| l IMF and then divided using Eq.[8] 
This approach, however, does not reproduce the observed stel- 
lar IMF as can be seen in Fig. [2] (thin lines), which displays the 
stellar IMFs produced in our simulations (filled squares, repre- 
senting simulation with initially 100% binaries). For compari- 
son, the standard IMF was included by setting up simulations 
using random pairing of the binary components instead of using 
Eq. [8] to pair the binary components. Random pairing has been 
shown to reproduce the stellar IMF ( |Leinert et al.|1993) and is 
here therefore used to represent the standard IMF. Using Eq. [8] 
to pair the binary components generates too few low-mass stars 
in the range 0.08M Q - O.16M and too many intermediate mass 
stars in the mass range 0.15 - 0.5M o , where the deviations be- 
come larger the higher the initial binary frequency is. The num- 
ber of stars generated with masses exceeding 0.5M G reproduces 
the [Kroupa (2002) IMF quite well, although their numbers are 
slightly under-predicted. 

The thick lines in Fig.[2]show the initial binary-system mass 
(t)hi„i.-M ■ j.M ) functions for the simulations used in this 

ny s-Jfip nm tiw sec ' 

study and additionally the initial binary-system mass function 
if the binary components are drawn randomly from the stellar 
IMF. The binary system mass function in both cases is very sim- 
ilar. We do not expect the differences in the single-star IMFs to 
affect the evolution of the binary populations because it depends 
on the system masses as we show in Sec. [4] 

The eccentricity of the binary systems in Sco OBII was 
found by Kouwenhoven et al. (2007 ) to be thermally distributed, 
in agreement with the eccentricities of long period binaries in 
the field. Following Kroupa & Burkert ( 2001[ ) who illustrated 
that the dynamical evolution is not capable of changing the ec- 
centricity distribution significantly, a thermal period distribution 
was applied to all of our simulations. 

Observations show that the binary frequency varies among 
low-mass stars in the field (35 



• 42%) ( [Fischer & Marcy|1992 
|Reid & G izis 1997), Sun-like stars in the solar neighbourhood 



100%) (Duche ne|1999| l. We account for the large spread among 
the different environments by investigating cluster models ini- 
tially with 30%, 50%, 75%, and 100% binary frequency. 

To improve the statistical significance of our results, for each 
initial binary frequency 100 simulations were performed. Only 
those that fulfilled the following conditions were considered in 
the subsequent diagnostic process: 

- after 1 Myr of cluster evolution, the most massive stellar sys- 
tem in the simulations has to be within lpc of the cluster 
centre. 

- the cluster density profile after 1 Myr of cluster evolution 
has to agree with the cluster density profiles observed by 
|Hillenbrand| ( [T997l l and |McCaughrean et al.| l |20Q2) . 

4. Results 

Observations of binary populations in different environments 
usually use the binary frequency as the main observable of 
a binary population since it can easily be compared be- 
tween different environments and different observing cam- 
paigns. Theoretical studies of the evolution of binary popula- 
tions are unaffected by from this shortcoming but usually also 
use this simple observable for reasons of comparison with ob- 



servational data (see for example Kroupa||1995a[ |Kroupa et al 



200T]|Kroupa & Bouvier 2003 Ivanova et al.|2005||Hurley et al. 



|2007i|Sollima|2008) . If we also do this here, we obtain the result 
described by Fig.|3£. It shows the evolution of the normalized bi- 
nary frequency during the first three Myr of cluster evolution in 
the simulations. The normalization is performed by dividing the 
given binary frequencies by the initial binary frequency and al- 
lows direct comparison of the evolution of the binary frequencies 
for simulations with different initial binary frequencies. 

As expected from previous studies ( ] Krou pa 1 995 aj [Kroupa] 
|& Bou vier 200"3] |Parker et al.||2009"] ), dynamical evolution re- 
duces the binary frequency with time in all simulations. This can 
be explained by the destruction of wide binary systems by means 
of three body encounters (see |Kroup a 1995a b , for more details). 
Apart from this general trend, Fig. [3^ shows that the evolution of 
the binary population depends on the initial binary frequency in 
the sense that the higher the initial binary frequency, the higher 
the fraction of binaries that are destroyed during cluster evolu- 
tion. This finding could lead to the conclusion that the evolu- 
tion of binary populations is stronger in clusters with a higher 
initial binary frequency. An analogous dependence was found 



(~ 60%) ( |Duquennoy & Mayor 1991 1, and the high binary 



frequency in low-mass star-forming regions such as Taurus (- 



by |Sollima (2008 1, who investigated the evolution of the binary 
population in globular star clusters using an analytical approach. 

In short, the evolution of the binary population in star clus- 
ters seems to depend on the initial binary frequency. However, 
because the calculation of the binary frequency depends on the 
number of single stars (see Eq. [2 1 it is probably not an unbiased 
measure of the evolution of the binary population. We thus in- 
troduce the normalized number of binary systems Nb (see Eq.[T]i 
as the proper quantity and show its evolution in Fig.[3[3. 

As expected from Fig. [3^, the normalized number of binaries 
also decreases with time. However, in contrast to the normalized 
binary frequency the evolution of the normalized number of bi- 
nary systems with time does not depend on the initial binary fre- 
quency. The difference among the data of only 1 % is of the order 
of the statistical error in the simulation. How can this difference 
in the evolution of the binary frequency and the binary popula- 
tion be explained? It is the evolution of the single star population 
in the star clusters that significantly affects the evolution of the 
binary frequency. Its evolution is show in Fig. [3J:. 
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Fig. 3. Evolution of the normalized binary frequency, normal- 
ized number of binary systems and normalized number of single 
stars with time. Crosses, stars, open, and filled squares denote 
simulations with initially 30%, 50%, 75%, and 100% binaries, 
respectively. 



In contrast to the binary frequency and the number of bi- 
naries, the number of single stars in the star cluster obviously 
increases with time because the dissolution of a binary system 
creates two new single stars. The evolution of the normalized 
number of single stars, however, depends even more strongly on 
the initial binary frequency than the evolution of the binary fre- 
quency. 

We show by means of simple arithmetic why the evolution of 
the single star population depends on the initial binary frequency 
while the evolution of the binary population does not. First, we 
assume that the evolution of the number of binary systems can 
be described by a function a(t,bf, . . . ) that may depend on the 
time of evolution, initial binary frequency, and other parameters, 
so that Nb(t) - a(t,bf, . . . )N b (0), where Nb(t) is the number of 
binary systems at time t. Dividing this by the initial number of 



binary systems we obtain the evolution of the normalized num- 
ber of binary systems 



N b {t) = 



a(t,b f ,...)N b <SS) 
N b (0) 



= a(t, b 



(9) 



If we also assume, based on the findings of Fig. [3]3, that the evo- 
lution of the normalized number of binary systems does not de- 
pend on the initial binary frequency, the function a can be ex- 
pressed as a(t,bf, . . . ) ~ a(t). Because the dissolution of a bi- 
nary system leads to the generation of two new single stars, the 
evolution of the number of single stars can be written as 



N s (t) = N.M + 2 [l - a(t,b f , ...)] N b (0), 



(10) 



where the first summand represents the initial number of sin- 
gle stars and the second summand gives the number of newly 

bf(P) 



formed single stars. Applying the relation Nb(0) = 



-N s (0) 



i-bf(0) 

and dividing by the initial number of single stars N s (0), we find 
that the evolution of the normalized number of binary systems is 
given by 

b f (0) 



N,(t) = 1 + 2(1 - o(0) 



1-bfQO) 



(ID 



Hence, in contrast to the normalized number of binary systems 
(Eq.[9]l, which depends only on a, the normalized number of sin- 
gle stars depends explicitly on the initial binary frequency in the 
cluster. Consequently, using Eqs.|9]and 11 we derive an analo- 
gous dependence for the normalized binary frequency 



b f (f) 



lb 



N b (t) 



1 



: '" N s (t)+N b (t)b f (0) 

a(t) 1 

b f (0) v ' 



(12) 



In summary, our simple model (Eqs. 9p2 1 describes the evo- 
lution presented in Fig. [3] very well: it is the evolution of the 
single star population that is very sensitive to the initial binary 
frequency, while the binary population evolves independently of 
the initial binary frequency. The normalized number of binary 
systems (Eq. [TJ introduced here is therefore the most simple 
measure of the evolution of a binary population independent of 
the initial binary frequency. 



As a consequence, we can use Eq. 12 to calculate the initial 
binary frequency in young ONC-like star clusters from only a 
few observables. Solving Eq.fTJlfor bj(Q) reveals that 



M0) = 



bf(t) - a(t) + 1 
ait) ' 



(13) 



This finding has a fundamental consequence: if we are able 
to determine the explicit form of the function a(t,bf, . . . ), we 
are able to predict the past and future number of binary sys- 
tems from the present number of binaries in any cluster environ- 
ment. If we again assume that a{t,bj, . . . ) — a{t), we find that 
ttftt(0 = at b +ct+\ with a = -0.35 + 0.04, b = 0.89 ± 0.01, and 
c - 0.27 + 0.04, provides an excellent fit to the data in Fig. [3}?. 
However, we emphasize that our assumptions neglects the po- 
tential effects of other parameters such as the stellar density or 
distribution of binding energy, which will be investigated in a 
forthcoming paper. 

Nonetheless, for illustrative purposes we assume that the 
ONC has cluster age of lMyr |Hillenbran"d| \\991) and obtain 
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a(lMyr) = 0.92 using the given fit. Assuming a present-day 
binary frequency of 60%, motivated by observations that find a 
consistency between the binary populations of field G dwarfs 
(|Duquennoy & Mayor|[l99"T) and the ONC ( |Petr et al.|[T998| 
|Kohler et al.||2006[ ), we derive an initial binary frequency of 
74%. This is clearly below the binary frequency observed in 
young sparse star-forming regions as Taurus -Auriga ( |Leinert| 
|et al.|1993||Ghez et al.|19 93) and suggests that the initial binary 
frequency depends on the star-forming environment as stated by 
|Durisen & Sterzik| ( |T994) . 

We note that this result is based on our specific cluster model, 
which does not include all the physical processes present in 
real star clusters. Thus, further investigations may alter the de- 
termined initial binary frequency to some extent. For example, 
including the gas the young clusters are embedded in would 
probably require a higher initial density since star clusters ex- 
pand significantly after gas expulsion (e.g. Baumgardt & Kroupa 
|2007| >. Higher densities should increase the fraction of dynami- 
cally destroyed binary systems during cluster evolution, result- 
ing in larger values of a. This means that the ONC could have 
been born with a higher initial binary frequency than estimated 



from our current simulations (as demonstrated by Kroupa et al. 
[2001) . 

We next use our homogeneous data set to investigate why 
the binary population evolves independently of the initial bi- 
nary frequency. Figure [4] shows the fraction of binary systems 
per semi-major axis bin that is destroyed after lMyr of cluster 
evolution. It is evident that only binaries with semi-major axis 
wider than about m 30AU are destroyed during cluster evolu- 
tion. This dependence on the semi-major axis was previously 
reported in other studies of the evolution of binary populations 
(e.g. Kroupa 1995b) and can be explained by the outcome of 
three body encounters (for details see Heggie 1975a i: to dissolve 



a binary system in a three body interaction, the kinetic energy of 
the intruding star must exceed the binding energy of the binary 
system, or, in other words, the velocity of the intruder must be 
larger than the orbital velocity of the binary system. With this, 
we estimate a lower limit to the semi-major axis of destroyed 
binaries in the cluster. Following Kroupa & Bouvier ( 2003), we 
estimate the "thermal period" log 10 P± « 4.75 at which destruc- 
tion still is possible from 



log 10 P th = 6.986 + log 10 m sys - 3 log 10 cr 3D 



(14) 



using the least massive equal-mass system m sys = O.16M and 
assuming a typical relative velocity equal to the cluster veloc- 
ity dispersion in our simulations cr 3£) ss 3km s . Our result, 
log 10 P t h ~ 4.75d, corresponds to a minimum semi-major axis 
of about 53AU. 

This estimate is in good agreement with the results of our 
simulations in Fig. [4] However, in our simulations, there is no 
sharp transition between stable and destroyed systems because 
the thermal period P,i, depends strongly on the system mass 
(Eq. [T4"| ). Since the periods of the binaries in the simulations 
have been sampled independently of the binary system masses 
the fraction of binaries destroyed by the dynamical evolution of 
the cluster decreases with binary mass. This is demonstrated in 
Fig. [5] The fraction of destroyed binaries clearly decreases with 
binary mass. 

The solid line in Fig. [5] represents a fit of the function 
fQn sys ) = alog 10 OT sys + b to the data of the simulations with 
initially 75% binaries in the range from ;« sys = 0.1 6M up to 
1OM . For higher masses, the results suffer from low number 
statistics and have therefore been excluded from the fit. We note 




1 10 100 1000 10000 

semi-major axis [AU] 

Fig. 4. Fraction of binary systems per semi-major axis bin that is 
destroyed after lMyr of cluster evolution. Crosses, stars, open, 
and filled squares denote simulations with initially 30%, 50%, 
75%, and 100% binaries, respectively. 
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Fig. 5. Fraction of binaries per binary system mass bin that are 
destroyed after lMyr of cluster evolution. The solid line is the 
fit of the function f{m sys ) - alog I0 ra sys + b to the data with 
initially 75% binaries. Crosses, stars, open, and filled squares 
denote simulations with initially 30%, 50%, 75%, and 100% bi- 
naries, respectively. 



that f(m sys ) also fits the data of the other simulations demonstrat- 
ing that the evolution of the binary population is independent of 
the initial binary frequency. 

The strong mass dependence of the binary destruction is fur- 
ther illustrated in Fig. [6^. It shows the evolution of the normal- 
ized number of binaries with primary masses either exceeding 
2M Q (thin lines) and lower than 2M (thick lines) in clusters 
with various initial binary frequencies. Clearly, a larger fraction 
of low-mass binaries is destroyed during the dynamical evolu- 
tion of the cluster. 

The boundary of 2M Q between high and low-mass primaries 
was chosen in accordance with the study of Kohler et al. (2006), 
who performed a high spatial resolution infrared study of the 
binary population in the ONC. They found the frequency of 
binaries with high-mass primaries to be significantly higher. 
Figure [6^ implies that this finding can be explained by the dy- 
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namical evolution of the cluster and does not depend on the ini- 
tial binary frequency. 

However, a consistent comparison of our simulation and the 
observations by Kohler et al. ( 2006 ) must focus on the frequency 
and not the number of binaries. Figure^ shows the evolution of 
the frequency of binaries with high and low-mass primaries us- 
ing the same symbols as in Fig. [6^. Clearly the binary frequen- 
cies of the high and low-mass primaries evolve quite differently. 
The former decreases much more slowly with time than the latter 
and does not depend on the initial binary frequency. In contrast, 
the evolution of binaries with low-mass primaries does depend 
on the initial binary frequency. 

The differences in the evolution are explained by two mecha- 
nisms. First, binaries with low-mass primaries are preferentially 
destroyed by dynamical interactions in the cluster because of 
their lower binding energy (Fig. [6^). Secondly, the destruction 
of a binary by a high-mass primary can result in the generation 
of two new single high-mass stars or a high and a low-mass star. 
The latter possibility means that the number of single low-mass 
stars increases, hence the frequency of binaries with low-mass 
primaries decreases relative to the former case. Simultaneously, 
the frequency of binaries with a high-mass primary does not de- 
crease as fast as implied by Eq. 12 because of the smaller in- 



crease in the number of single high-mass stars. Hence, the pos- 
sible generation of a high- and a low-mass star in the destruction 
of a binary with a high-mass primary accelerates the reduction 
of the binary frequency of low-mass stars and decelerates the 
reduction of the binary frequency of high-mass stars. 

We emphasize the different evolution of binaries with high 
and low-mass primaries in Fig.[6j;. Here we show the difference 
between the binary frequencies of high and low mass stars over 
time, Ahibf(t), relative to its initial value Ahi(f)- This normaliza- 
tion accounts for the possible initial difference of the binary fre- 
quencies of high and low-mass stars, so that all values are zero 
initially. 

In all simulations, the difference of the binary frequency of 
high and low-mass stars increases with time because the binary 
frequency of low-mass stars decreases faster than the binary fre- 
quency of high-mass stars (Fig. ^p). However, the amount of 
increase depends strongly on the initial binary frequency in the 
cluster with the lowest initial binary frequency (30%) resulting 
in the smallest increase of less than 2% after 3Myr while in 
simulations with the maximum initial binary frequency (100%) 
the difference increases by more than 10%. Under the assump- 
tion that initially the frequency of binaries with low- and high- 
mass stars was the same in the ONC, the observational data from 
Kohler et al. (2006) imply that the binary frequency in the ONC 
initially was very high with values well above 75%. If the initial 
binary frequencies were lower initially the observed difference 
must have been part of the initial conditions of the binary popu- 
lation. 

5. Conclusions 

We have studied the dynamical evolution of the binary popu- 
lation in young dense ONC-like star-clusters using two differ- 
ent measures: the (normalized) binary frequency (Eq.|2]) and the 
normalized number of binary systems (Eq. [TJ. As expected, the 
evolution of the binary frequency depends heavily on the ini- 
tial binary frequency: the higher the initial binary frequency in a 
cluster, the stronger the decrease in the binary frequency during 
cluster evolution. However, our main result is that the evolution 
of the normalized number of binaries is independent of the ini- 
tial binary frequency. During the first Myr of cluster evolution, 
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Fig. 6. Panel a: Evolution of the normalized number of binary 
systems with primary masses exceeding 2M Q (thin lines) and 
binaries with masses below 2Mo (thick lines) with time. Panel 
b: Evolution of the normalized binary frequencies of binaries 
with high- (thin lines) and low-mass (thick lines) primaries with 
time. Panel c: Evolution of the difference of the binary frequency 
of high- and low-mass stars with time. The difference at each 
time-step has been reduced by the initial difference of the binary 
frequencies to improve the readability. Crosses, stars, open, and 
filled squares denote simulations with initially 30%, 50%, 75%, 
and 100% binaries, respectively. 



the number of binary systems decreases by about 6 - 9% inde- 
pendent of the initial binary frequency in the clusters. 

The reason for this discrepancy between the two measures 
is the dependence of the binary frequency on the population of 
single stars. Since the destruction of a binary system leads to 
the generation of two new single stars, the single star population 
can easily dominate the evolution of the binary frequency in a 
cluster, especially in the case of high initial binary frequencies 
(Eq. 12 1. Hence, the number of binaries - unbiased by the evo- 
lution of the single star population - represents the proper mea- 



sure of the evolution of the binary population. We can therefore 
conclude that the evolution of the binary population in ONC- 
like stars clusters is independent of the initial binary frequency, 
which is in striking contrast to what would be inferred from the 
evolution of the binary frequency. 

This independence of the number of binaries has allowed us 
to determine the initial binary frequency and predict the evo- 
lution of the binary population in young ONC-like stars clus- 
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ters from present-day observational data. For example, taking 
the observed binary frequency of the ONC of ~ 60% and its 
age of lMyr ( |Petr et al.|1998||Hillenbrand|1997] l, and estimating 
Q'(lMyr) = 0.92 from our fit to the simulations, Eq. [T~3| provides 
an initial binary frequency of 74%. This binary frequency is sig- 
nificantly lower than the binary frequency observed in sparse 
star-forming regions such as Taurus-Auriga (e.g. |Leinert et al.| 
1993 Ghez et aTf l993 ) and implies that the initial binary fre- 
quency depends on the star-forming region. This finding agrees 
well with the results of |Durisen & Ster zik ( 1994), who claimed 
that the primordial binary frequency in star-forming regions de- 
pends on the temperature of the natal molecular cloud in which 
the stars are forming of. However, different initial conditions of 
the clusters (e.g. higher initial stellar densities) could possibly 
increase the amount of destroyed binary systems, which would 



further increase the estimated initial binary frequency (Kroupa 
et al.||200T] >. It is not trivial to quantify these effects, so further 
numerical studies are required. 

In addition, our simulations show that the evolution of the bi- 
nary populations depends crucially on the masses of the primary 
stars: the lower the mass of the primary, the more likely that the 
binaries are destroyed. Over time, this leads to a higher binary 
frequency for stars with masses exceeding 2M Q than for lower 
mass stars. This is even so in the case where the binary frequen- 
cies of high and low-mass stars were initially the same. This 



finding is in good agreement with the observations of |Kohler 



et al. 



( 2006 ), who found a significantly higher binary frequency 



for high-mass stars than for low-mass stars in the ONC. Our re- 
sults suggest that the difference in the binary frequencies of high- 
and low-mass stars are probably not the outcome of the binary 
formation process but rather the result of the dynamical evolu- 
tion with time. A possible test of this hypothesis would be to 
compare the maxima of the period/semi-major axis distributions 
of high- and low-mass stars in young dense star clusters. Since 
the probability of destroying a binary in a three body encounter 
depends on the binding energy and therewith on the separation 
and mass of the binary, the maximum of the period and semi- 
major-axis distributions for massive stars should be at longer 
periods and wider semi-major axis than for lower mass stars. 
However, this test is not easy to be earned out since observations 
of the binary population in young dense star clusters usually only 
cover a small range of periods/separations of the binaries, and it 
is often not even possible to determine the physical orbital pa- 
rameters based on the observations. 

The simulations presented here were performed for ONC- 
like clusters only and used the specific binary population model 
of Kouwenhoven et al. ( 2007 1. We aim to generalise the present 
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